*------------------------------------------------------------------------------*
			*** FIGURE S12b. HISTOGRAMS OF RESIDUALS WITH CONTROLS ***
*------------------------------------------------------------------------------*

use "${tempdata}kreise_complete_all.dta", clear

*------------------------------------------------------------------------------*

***Residuals plot with controls
qui reg PM10_30km_10ug O3_30km_10ug temp_30km humid_rel_30km wind_speed_30km precip_30km i.kreis_code i.election_date, vce(cluster kreis_code), if samplePM10_allcontrols == 1
predict res_wc, residuals

*Histogram
histogram res_wc, fcolor(navy) lcolor(ltbluishgray%20) lwidth(vthin) graphregion(fcolor(white))
graph export "${outputs}figure_S12b.pdf", as(pdf) replace

*------------------------------------------------------------------------------*

clear

exit
